A weak comparison principle for 
reaction-diffusion systems 



Jose Valero 

Centro de Investigation Operativa 

Universidad Miguel Hernandez, Avda. Universidad, s/n. - Elche - 03202 
e-mail: jvalero@umh.es 



Abstract 

In this paper we prove a weak comparison principle for a reaction-diffusion system 
without uniqueness of solutions. We apply the abstract results to the Lotka-Volterra 
system with diffusion, a generalized logistic equation and to a model of fractional-order 
chemical autocatalysis with decay. Morever, in the case of the Lotka-Volterra system 
a weak maximum principle is given, and a suitable estimate in the space of essentially 
bounded functions L°° is proved for at least one solution of the problem. 
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1 Introduction 

Comparison results for parabolic equations and ordinary differential equations are well known 
in the literature. One of the important applications of such kind of results is the theory of 
monotone dynamical systems, which leads to a more precise characterization of w-limit sets 
and attractors. In the last years several authors have been working in this direction (see, 
for example, [3], [TO], [TO], [20], [21] for the deterministic case, and [I], [2], [8], [TOJ for 
the stochastic case). In all these papers it is considered the classical situation where the 
initial-value problem possesses a unique solution. 

However, the situation is more complicated when we consider a differential equation for 
which uniqueness of the Cauchy problem fails (or just it is not known to hold). Let us 
consider an abstract parabolic problem 

-£ = A(t,u(t)), r<t<T } (1) 
u (r) = u T , 

for which we can prove that for every initial data in the phase space X (with a partial order 
<) there exists at least one solution. 
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If we try to compare solutions of ([I]) for two ordered initial data u\ < v%, then we can 
consider a strong comparison principle and a weak one. 

The strong version would imply the existence of a solution u\ with u\ (r) = u\ such that 

ui{t) <u 2 (t) for t E [r,T], (2) 

for any solution u 2 with u 2 (r) = u 2 , and, viceversa, the existence of a solution u 2 with 
u 2 (r) = such that (j2j) is satisfied for any solution U\ with n x (r) = it*. This kind of result 
is established in [S] for a delayed ordinary differential equations, defining then a multivalued 
order-preserving dynamical system. 

The weak version of the comparison principle says that if u\ < u 2 , then there exist two 
solutions ui,u 2 of ([1]) such that u\ (r) = u\, u 2 (r) = ul, and ([2]) holds. 

There is in fact an intermediate version of the comparison principle, which says that if 
we fix a solution ui of ([T]) with u\ (r) = u\, then there exists a solution u 2 with u 2 (r) = v? T 
such that (T5]) is satisfied (and viceversa). This is proved in [6] for a differential inclusion 
generated by a sub differential map. 

In this paper we establish a weak comparison principle for a reaction-diffusion system 
in which the nonlinear term satisfies suitable dissipative and growth conditions, ensuring 
existence of solutions but not uniqueness. This principle is applied to several well known 
models in Physics and Biology. Namely, a weak comparison of solutions is proved for the 
Lotka-Volterra system, the generalized logistic equation and for a model of fractional-order 
chemical autocatalysis with decay. Morever, in the case of the Lotka-Volterra system a weak 
maximum principle is given, and a suitable estimate in the space of essentially bounded 
functions L°° is proved for at least one solution of the problem. 

We note that in the papers [TT] , [T2] the existence of a global attractor is proved for such 
kind of react ion- diffusion systems. In a near future we will apply these results to obtain 
theorems concerning the structure of the global attractor. 

2 Comparison results for reaction-diffusion systems 

We shall denote by |-| and (•,•) the norm and scalar product in the space R m , m > 1. Let 
d > be an integer and Q C M N be a bounded open subset with smooth boundary. Consider 
the problem 

du 

— - aAu + f(t, u) = h(t, x), (t, x) e (r, T) x Q, 

uUan = 0, ( 3 ) 
u\ t=T = u T (x), 



where r, T G R, T > r, x E Q, u = (u 1 (t, x) , u d (t, x)) , / = (f 1 , f d ) , a is a real 

d x d matrix with a positive symmetric part — — 01, > 0, h E L 2 (t,T; {I? (tt)) d ). 
Moreover, / = u), f d (t, u)) is jointly continuous on [r,T] x R d and satisfies the 
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following conditions: 

d d 

J2\f%u)\^ KC^l + ^T), (4) 



i=l i=l 



(f(t,u),u)>aJ2\uT-C2, (5) 

where Pi > 2, a, Ci, C 2 > 0. 

Let H = (L 2 (tt)) d , V = (H% (tt)) d , and let V be the dual space of V. By ||-|| , we 
denote the norm in H and V, respectively. For p = ...,p<z) we define the spaces 

L p (Q) = L P1 (fi) x • • ■ x L Pd (fi) , 
L p (r, T; L p (ft)) = L Pl (r, T; L Pl (ft)) x • ■ ■ x L Pd (r, T; L Pd (ft)) . 

We take q = (q u q d ), where i + i = 1. 

We say that the function u (-/is a weak solution of (ED if u G L p (r, T; L p (ft))n L 2 (r, T; V)n 

C([r, T]; ), — G L 2 (t, T; V) + (r, T; L« (ft)) , u (r) = u T , and 

y (^,t\dt + J J (V (aw) ,V£)dxdt + J J (f (t,u) ,£) dxdt = J J (h,£)dxdt, 

(6) 

for all £ G L p (r, T; i7(ft))n L 2 (r, T; V), where (•,•) denotes pairing in the space V + L q (ft), 
and (Vw, Vv) = Eti (Vtt*, W) . 

Under conditions (jljl-flHJ) it is known [7J p. 284] that for any u T £ H there exists at least 
one weak solution u = u(t,x) of (ED, and also that the function £ !->■ ||w(t)|| 2 is absolutely 

continuous on [r,T] and — ||m(;£)|| 2 = 2 / for a. a. t G (t,T). 

dvu 

Denote r = (r 1; r^) , = max{l; iV (l/<& — 1/2)}. Any weak solution satisfies — G 
L 9 (r,T;#- r (ft)) and 

L q {0,T;H~ r {Q)) = L qi {0,T;H- ri {Q)) x ••■ x L qd {0,T; H~ rd {Q)). 

If, additionally, we assume that that / (t, w) is continuously differentiable with respect to 
u for any t G [r, T] , w G M d , and 

(/«(*, u)w, w) > -C 3 (t) \w\ 2 , for all w, u G M d , (7) 

where C^-) G L 1 (r,T), C3 (£) > 0, the weak solution of d2J) is unique. Here, /„ denotes the 
jacobian matrix of /. 

We consider also the following assumption: there exists R > such that 

f(t,u)>f(t,v), (8) 
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for any t G [r, T] and any u, v G M. d such that u l = v % and < f- 7 if j ^ i, and \u\ , \v\ < Rq, 
which means that the systems is cooperative in the ball with radius Rq centered at 0. 
Consider two problems 



du 

— -oAu + /i(t,u) = h(t,x), (t,x) G (r,T) x Q, 

u\xedn = 
u\t= T = u T (x) 



(9) 



(9u 

— - aAu + f 2 (t, u) = h 2 (t, x), (t, x) G (r, T) x J], 

■u|ccesn = 0, 

U\ t =r = U T (x), 



(10) 



where fj are jointly continuous on [r, T] x ~R d . Among conditions (jlj)-© and ©-([8]) we shall 
consider the following: 

hi (t, x) < hi (t, x) , for a. a. (t, x) , (11) 
fx (t, u) > f 2 (t, u) , for all t, u. 

Lemma 1 If fj satisfy (0)-|3J) and 07]) . then the constants Pi has to be the same for /i and 

h- 

Proof. Denote by p\, oP , C{, C 3 2 the constants corresponding to fj in (JIJ-©. By contra- 
diction let, for example, p\ > p\. Take the sequence u n = {u\, 0, 0) , where u\ — > +oo as 
n — )• oo. Then by (BJ-fEJ), ( 031) and Young's inequality we have 



a 2 \v} n \ Pl - C\ < (f 2 (t,u n ),u n ) < (fi(t,u n ),u n ) 



p}_-i 



= fl(t, Un )u 1 n < (cfti + Kl 1 *)) Pl <<C(1 + | W ^). 

But > implies the existence of n such that a 2 — Cf > C(l + |m^| p i), which is a 
contradiction. Hence, < p\. 

Conversely, let p\ < p\. Then we take u n = (m^,0, ...,0) with v} n — > — oo as n — ^ oo, so 
that 

a 1 \u l n \ Pl - C\ < (fi(t,u n ),u n ) < (f 2 (t,u n ),u n ) 



= f\ (t,u n )u l n < (C?(l + |t£|rf))~*~t*? < C(l + 

As before, we obtain a contradiction, so = p\. 

Repeating similar arguments for the other p\ we obtain that p\ = p 2 for i = 1, d. 
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We recall [TT] that under conditions (jlj)-© any solution u (■) of (jUJ) satisfies the inequality 

t d t t 

\u(t)\\ 2 +2/3 J \\Vu(T)\\ 2 dT+aJ2 J \\u\r)\\% m dr < \\u(s)\\ 2 +C J QMr^+^dr, (12) 



8=1 



for some constant C > 0. Of course, the same is valid for any solution of fllOp . From (1121 , 
for any T > r we obtain 

T 

||«(*)|| 2 < \K\\ 2 + C J (Wh^r)]] 2 + l)dr = K 2 (\\u T \\,r,T) for all r < t < T. (13) 

T 

We shall denote by ui(t) the solution of ([9]) corresponding to the initial data u l Tl and by 
u 2 (t) the solution of (flOl) corresponding to the initial data u 2 . Also, we take u + = max{«, 0} 
for uGi 

We obtain the following comparison result. 

Theorem 2 Assume that fj,hj satisfy OU-fGP, (?]) and (TJIj). // < u 2 and we suppose 
that f 2 satisfies (TJ) with R 2 > 2 m&x{K 2 (\\u% r, T) , if 2 (||u 2 ||, r, T)}, where K r, T) 

zs ta/cen from < f73|) . we /iai>e Wi (t) < «2 (0; / or a ^ * ^ [ r > 

Remark 3 The results remains valid if, instead, f\ satisfies (Gj) with Rq > K 2 (\\u^.\\,T) . 

Remark 4 If f 2 satisfies (Gj) for an arbitrary Rq > (that is, in the whole space M. d ), then 
the result is true for any initial data u\ < u 2 . 

Proof. Let g 2 (t, u) = f 2 (t, u) + C 3 (t) u. The function g 2 (t, •) satisfies ([7]) with C 3 = 0. For 

any Ui,u 2 G M d define v 2 (ui,u 2 ) by 

v i = (u\,Hu[>u 2 , 
2 \u\, if u\ <u l 2 . y ' 

Note that U\ — v 2 = {u\ — u 2 ) + and v 2 — u 2 = — (u 2 — «i) + , so v 2 < u\ for all %. For the 
function (u\ — u 2 ) + we can obtain by (fTTj) and the Mean Value Theorem that 



2dt 



(m - u 2 ) + \\ 2 + a - u 2 ) + \\ 2 v < - / (/i (t,ui) - f 2 (t,u 2 ) , (m - u 2 ) + ) dx 

Jn 

< ~ I (/a Mi) - / 2 (*, u 2 ) , (ui - u 2 ) + ) dx- (f 2 (t, v 2 ) - f 2 (t, u 2 ) , {u x - u 2 ) + ) dx 

Jn 

{g2u {t,v(t,X,U 1 ,U 2 }) {U X - V 2 ) , (u X - U 2 ) + ) dx + C 3 (t) || (li! - M 2 ) + || 2 

(/2 U2) - /a U2) , (wi - m 2 ) + ) dx 
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-/ (g2u(t,v(t,x,u 1 ,u 2 ))(ui-u 2 ) ,{u 1 -u 2 ) ) dx + C 3 (t) II («i - « a ) ||~ 

(/2 (*, ^2) - /2 (*, W2) , ("i - M 2) + ) <fe 

<C 3 (t)||K-n 2 ) + || 2 - / (f 2 (t,v 2 )-f 2 (t,u 2 ),(u 1 -u 2 ) + )dx. 

Jn 

For all t we have 

(/ 2 (t, v 2 ) - f 2 (t, u 2 ) , ( Ul - u 2 ) + ) = (fi (t, V2) - tt (t, u 2 )) (u\ - 4) + , (15) 

where u\ — u l 2 > 0, for i G J, and u\ — u l 2 < if i J. For any i G J we have that u| = w|, 
and then by v 2 < u 2 , \v 2 \ 2 < |mi| 2 + \u 2 \ 2 , (IT51) and (jSj) we get 

n(t,v 2 )-n(t,u 2 ) >o. 

By Gronwall's lemma we get 



|K (t) -n 2 (t)) + || 2 < {ul-ul) 



' e f*2C 3 (s)ds = q 



Thus |(wi — ^2) + || = 0, which means that u\ (x,t) — m 2 (x, t) < 0, for a. a. x G f2, and all 
ie{l,..,d},te [r,T]. ■ 

Remark 5 In the scalar case, that is, d = 1, condition (TJJ) trivially satisfied. 

When condition (JTj) fails to be true, we will obtain a weak comparison principle. 
Define a sequence of smooth functions ip k : R+ — > [0, 1] satisfying 

1, if < s < fc, 

^ (s) = { < V fe (s) < 1, if A; < s < k + 1, (16) 
0, if s > k + 1. 

For every fc > 1 we put /^(t, u) = ^ (M) / l (i, u) + (1 — V'fc (M)) ( M ) ; where (7* (it) = 
|it i | IH ~ 2 ii i . Then / fe G C([r,T] x R d ; R d ) and for any A > 0, 

sup sup |/fc(t, m) — fit, u) I —j- 0, as k — >■ 00. 

te[r,T] |«|<A 

Let p £ : R d -»■ R+ be a mollifier, that is, p t G Cg° (R d ; R) , supp p e C B e = {x E R d : 
\x\ < e}, L d p e (s) (is = 1 and p e (s) > for all s G R d . We define the functions 

/fc(^ M )=/ p t (s)f k (t,u- s)ds. 



Since for any k > 1 fk is uniformly continuous on [0, T] x [— k — 1, fc+1], there exist G (0, 1) 
such that for all u satisfying \u\ < fc, and for all s for which |it — s\ < e k we have 

sup \f k (t,u)-f k (t,s)\<\. (17) 

*e[r,T] fc 

We put f k (t,u) = f k k (t,u). Then f k (t, ■) G C°°(R d ; R" 1 ), for all * G [r, T], jfe > 1 . 
For further arguments we need the following technical result [T21 Lemma 2]. 

Lemma 6 Let f satisfy np~(EP. For all k > 1 the following statements hold: 

sup sup \f k (t,u) - f{t,u)\ ->■ 0, as jfe ->■ cx), V A > 0, (18) 

te[r,T] |u|<A 
c? d 

J2 \f kl (t,u)\^ <D 1 (l + 53 (/*(«,«),«) >7X)l« < l ft - ^2, (19) 

i=l i=l i=l 

(/ u fc (t,u) W , W ) > -D 3 (A;) |w| 2 , Vu.w, (20) 

where D 3 (k) is a non-negative number, and the positive constants Di, D 2 > C 2 , 7 do not 
depend on k. 



Consider first the scalar case. 

Theorem 7 Let d = 1. Assume that fj,hj satisfy (^-(EJ) and < f77)j . < w 2 , £/iere exist 

two solutions u±,U2 (of ^ and (TQ), respectively) such that u\ (t) < u 2 (t), for all t G [r, T]. 

Proof. For the functions fj we take the approximations f k (defined in Lemma [6]), which 
satisfy (|T8|) - (!20|) . and consider the problems 

— - aAu + f k (t, u) = hj (t, x) , (t, x) G (r, T) x Q, 

■\xEdCi = 0, 

u\ t=T = U T , 

for j = 1,2. Problem ( l2Ti) has a unique solution for any initial data w T G if. In view of 
Lemma [1] the constant p is the same for f 1 and / 2 . We note that 

hk {t,u)=i> k (H) /i (*, «) + (i - ^ k (\ u \)) \ u r 2 u 

> iJ k (\u\) h (t, u) + (1 - ^ (M)) \u\ p - 2 u = f 2k (t, u) . 

Then it is clear that f k (t, u) > f 2 (t, u) for every (t, u) . 

By Theorem [2] we know that as u\ < we have u\ it) < u\ (t), for all t G [r, T], for the 
corresponding solutions of (l2~Pj) . 
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In view of Lemma [6] one can obtain in a standard way that ffl2|) is satisfied for the 
solutions of (T2T1) with a constant C not depending on k and replacing a by 7. Hence, the 
sequences u) (■) are bounded in L°° (r, T; H) f)L 2 (r, T; V) nL p (r, T; L p (fi)). It follows from 

du k 



T9|) that /,&(-, (•)) are bounded in L q (r, T; L q (fl)) and also that {— — (•)} is bounded in 



L q (r,T; H r (fi)), where = max{l; (| — — )iV}. By the Compactness Lemma [H] we have 



that for some functions Uj = %(£, x), j — 1, 2: 

uj -> weakly star in L°°(r, T; if), (22) 

in L 2 (r,T;#), u){t) —> Uj(t) in if for a. a. t G (r, T), (23) 
u*(i,a;) -» Uj(t,x) for a.a. (t,x) G (r, T) x Q, (24) 
uj -> Uj weakly in L 2 (r, T; V), (25) 

-> weakly in L*(T,T;fr p (n)), (26) 

uj ->■ u 3 - weakly in L p (r, T; L p (fi)). (27) 
Also, arguing as in [T7J p. 3037] we obtain 

-> weakly in H for all t G [r, T]. (28) 

Moreover, by (JTHD and (j21J) we have f jk (t,u k (t,x)) ->■ fj(t,u(t,x)) for a.a. (i, a;) G (r,T) x 
and then the boundedness of fjky,u^(-)) in L q (r,T; L q (Q)) implies that fjk(-,Uj (■)) 
converges to /(-,«(•)) weakly in L q (r,T] L q (Q)) [14] . It follows that «i(-), W2( - ) are weak 
solutions of (jUJ) and ffTOl . respectively, with «i(r) = m*, u 2 (r) = u 2 . 
Moreover, one can prove that 

u){t) ->■ strongly in # for all i G [r, T]. (29) 

i 

Indeed, we define the functions Jjk{t) = \\u k (t)\\ 2 - C J{\\hj{s)\\ 2 + l)ds, Jj(t) = \\uj(t)\\ 2 - 

T 

t 

C J( \\hj(s) || 2 + l)ds, which are non-increasing in view of (1121) . Also, from (I2"3"j) we have 

T 

Jjkif) — > J j{t) for a.a. t G ij,T). Then one can prove that limsup^^ Jjfc(t) < J jit) for all 
i G [r,T] (see [11], p. 623] for the details). Hence, lim sup^^ < 11^(011- Together 

with (I2"gj) this implies (1291 (see again [TTJ p. 623] for more details). 
Hence, passing to the limit we obtain 

Mi it) < u 2 it) , for all t G [r,T]. 



Further, let us prove the general case for an arbitrary del 
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Theorem 8 Assume that fj,hj satisfy (0)-|3J) and Also, suppose that either f\ or f 2 

satisfies |2P for an arbitrary Rq > 0. Ifu\< u 2 , there exist two solutions Ui,u 2 (of and 
l[T(J\) . respectively) such that u\ (t) < u 2 (t), for all t 6 [t,T]. 

Proof. Let fi be the function which satisfies ([8]). We take the approximations f 2 
(defined in Lemma [6]), which satisfy (|T8|) - (!20|) . Then we consider problems fl2T|) . 
In view of Lemma [1] the constants pi are the same for f 1 and f 2 . We note that 

r lk (t, u) = ^ k (\ u \) fi (t, u) + (i - ^ k (|«d) it**!"- 2 

> ^ k (\ u \) f* (t, u) + (i - (H)) ju« p~ 2 = & a, u) . 

Then it is clear that /* (t, u) > f 2 (t, u) for every (t, u) . 

Using Lemma [6] it is standard to obtain estimate f|T3|) with a constant C not depending 
on k. Hence, the solutions u k (■) of (12~Tj) satisfy 

T 

114 CON 3 < IKf + Cj (IMrW + l)dr = K 2 (||<||,r,T) . 

T 

We note that 

/lib (t, u) = fi (t, u) , 

if \u\ < k, since in such a case ip k (\u\) = 1. Hence, if fc 2 > 2 max{K 2 ||, r, T) , i^ 2 (||n 2 1|, r, T)} 
the functions /i*. satisfy condition (jSJ) with i? = Therefore, for any £ 6 [r, T] and any 
m,d6 M d such that w l = v 1 and w- 7 < f- 7 if j 7^ z, and \u\ ,\v\ < k — 1, we have 

f*(t,u)= p €k {s)fi k {t,u- s)ds> p ek (s)fi k {t,v-s)ds = f$(t,u). 

Thus, if (k-lf > 2m&x{K 2 (\\u%T) , K 2 (||w 2 ||,T)}, the functions ft satisfy condition 
USD with R = k - 1. 

By Theorem [2] we know that as < w 2 , we have 14* (£) < i4(0> f° r a ^ ^ e [ r >^1> 
k > 1 + (2max{ir 2 (||^||, T) , K 2 (||w 2 ||, T)})s, for the corresponding solutions of f[2Tj) . 

Repeating the same proof of Theorem [7] we obtain that the sequences u\, u\ converge (up 
to a subsequence) in the sense of fl2"2"|) - fl2"9"|) to the solutions U\,u 2 of problem (j5J) and ffTUl) . 
respectively. Also, it holds 

Ux(t) <u 2 {t), for all t G [r,T]. 



In the applications we need to generalize this theorem to the case where the constant a 
can be negative. We shall do this when fi, f 2 have sublinear growth (that is, Pi = 2 for all 
i). Consider for ([3]) the following conditions: 

|/(t,u)|<Ci(l + |u|), (30) 
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(f(t,u),u) >a\u\ 2 -C 2 , (31) 

where «6l, and C±, C 2 > 0. 

Let fx, fa satisfy (!30|) - (l3T|) with constants a\ C\,C 3 2 , j = 1,2. Then if minja 1 ,*^ 2 } < 

0, we make in ([3]) the change of variable v = e~^u, where > — minja 1 , a 2 }. Hence, 

-pt 



(32) 



multiplying (J2J) and ( TlOl) by e 1 we have 
dv 

— - aAv + e~ pt fa{t, ePv) + 0v = e'^hxit, x), (t, x) G (r, T) x fi, 

v\ X £dn = 0, 
v\t= T = e~^ T u T (x), 

3d 

aAv + e-^fa (t, e^v) + 0v = e~^h 2 (t, x) , (t, x) e(r,T)xO, 



9* 



(33) 



e ^ r « T fx). 



It is easy to check that if v (t) is a weak solution of (l32l) . then it (t) = e^*i> (t) is a weak 
solution of ([9]) (and the same is true, of course, for ( 1331) and (TTU1) ). Conversely, if u(t) is a 
weak solution of (EJ), then t> (t) = e _/3 *u (t) is a weak solution of f l32l . 

The functions (t, v) = e~P l fj(t, e^v) + (3v satisfy (HJ)-® with Pi = 2 for all i. Indeed, 

\fj(t, v)\ << e-^Ci (1 + ef* \v\)+(3 \v\ < C{ (1 + \v\), (34) 

(fj(t, v),v)> e-^ifjit, e^v), e*v) + \v\ 2 (35) 
>{ot+0) \v\ 2 -Cl 

where a 3 + f3 > 0. 

Then we obtain the following. 

Theorem 9 Assume that fj,hj satisfy (3ty) - (3l\) and (Tl\) . Also, suppose that either fa or 
fa satisfies (EJ) for an arbitrary Rq > 0. If ' u\ < u 2 ., there exist two solutions Ux,u 2 (of (EJ) 
and / TTOj) . respectively) such that Ui (t) < u 2 (t), for all t G [r, T]. 

Proof. We consider problems ([32]Und ([33]). In view of (I54>(l3"5l) fa (t, v) = e'^fait, e^v) + 
f3v satisfy (@J-([5]). Also, defining hj (t,x) = e~^ t hj(t,x) it is clear that fTTTT) holds. Finally, 
if, for example, fa satisfies ([8]) for any Rq > 0, then it is obvious that for fa is true as well. 

Hence, by Theorem [8] there exist two solutions Vx,v 2 (of fl32]) and fl33l . respectively), with 
Vj (r) = e _/3r i4, such that v\ (t) < v 2 (t), for all t G [r, T]. Thus 

Ul (t) = e^vx (t) < e pt v 2 (t) = u 2 (t) , for t G [r, T], 

and Ui,u 2 are solutions (of ([9]) and ffTO]) . respectively, such that Uj (r) = u-?. ■ 

Remark 10 If fj satisfy {?]], then the solutions Ui,u 2 given in Theorem^ are unique for 
the corresponding initial data. 
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3 Comparison for positive solutions 

Denote = {u G M. d : u l > 0}. Let us consider the previous results in the case where the 
solutions have to be positive. Consider now the following conditions: 

The matrix a is diagonal, (36) 

h* (t, x) - f (t, u 1 , 0, u i+l , u d ) > 0, (37) 

for all i, a.e. (t,x) G (r, T) x fi and w- 7 > if j ^ i. Obviously, in the scalar case these 
conditions just mean that 

h(t,x)-f(t,0)>0, (38) 

for a.e. (t,x) G (r, T) x Q. 

It is well known (see [T2l Lemma 5] for a detailed proof) that if we assume conditions 
(J1J)-(|SJ) only for u G R+, and also (JTj) and (l36l) -( l3~Tj) . then for any u T > there exists a unique 
weak solution u (■) of (jHJ). Moreover, u (•) is such that u (t) > for all t G [r, T]. 

On the other hand, if we assume these conditions except (JTJ), then there exists at least 
one weak solution u (•) of (j3J) such that u (t) > for all t G [r, T] [12, Theorem 4]. Moreover, 
we can prove the following. 

Lemma 11 Assume conditions only for u G and also (3o))-(3l). Then 

there exists a weak solution it (•) of (TJ|), which is unique in the class of solutions satisfying 
u (t) > for all t G [r, T] . 

Proof. Let u\, U2 be two solutions with it, (r) = it r , z = 1,2, and such that -Uj (£) > for all 
t. Denote w (t) = Ui (t) — u 2 (t). Then in a standard way by the Mean Value Theorem we 
obtain 

X d f 

2dt "™ ® " 2 ~ ~ J Ul ^' ^ _ ^ ^' U2 ^' ^ ' ^ ^' ^ dX 
(/u v (t, %, Ui, u 2 )) w (t, x) , w (t, x)) dx 
<C 3 (t) \\w(t)f, 

where v (t, x, Ui, u 2 ) G L (ui (t, x) , u 2 (t, x)) = {aui (t, x) + (1 — a) u 2 (t, x) : a G [0, 1]}, so 
that, v (t,x,ui,u 2 ) > 0. The uniqueness follows from Gronwall's lemma ■ 

We prove also a result, which is similar to Lemma ED Denote by p\, a 3 , C{, C 2 the 
constants corresponding to fj in (JI|)-([5]). For problems (jHD and ( llOj) . respectively. Arguing 
as in the proof of Lemma [1] we obtain the following lemma. 

Lemma 12 If fj satisfy and < f77]j /or w G Et ; then p\ > pf for all i. 
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Theorem 13 Let fj,hj satisfy ([?]) and l[3b)) -(3l^. Assume that fj,hj satisfy 0-([3P and 
l[Tl\) for u G IRj. // < u\ < u 2 . and we suppose that f 2 satisfies (Ej) for u G with 
R 2 > 2max{K 2 (\\u 1 T \\,T,T),K 2 (\\u 2 T \\,r,T)}, where K r, T) zs taken from we 

have < Mi (t) < u 2 (t), for all t G [r, T], where U\ (■) , u 2 (•) are £/ie solutions corresponding 
to u\ and u 2 , respectively. 

Proof. As the solutions U\ (•), u 2 (•) corresponding to u\ and u 2 are non-negative, repeating 
exactly the same steps of the proof of Theorem [2] we obtain the desired result. ■ 

Remark 14 The results remains valid if, instead, f\ satisfies (EP with the same Rq. 

Remark 15 If f 2 satisfies (EP for an arbitrary Rq > (that is, in the whole space M. d ), then 
the result is true for any initial data < u\ < u 2 . 

We shall need also the following modification of Theorem [TBI 

Theorem 16 Let fj,hj satisfy ([?[) and l[3b)) -(3l^. Assume that fj,hj satisfy 0-([3P and 
HI]) for u G Let < u\ < u 2 . We suppose that f 2 satisfies 

ft(t,u)>f 2 (t,v)-e, (39) 

for any t G [r, T] and any u, v G such that u l = v l and u J < f J if j ^ i, and \u\ , \v\ < Rq 
with Rl > 2ma 1 x{K 2 (\\u 1 T \\,T } T),K 2 (\\u 2 T \\,T,T)}, where K {\\u{\\ , r , T) is taken from (E 
Then there exists a constant C (r, T) such that 

| |(M*) ~M*)) + || < C(r,T)e, for all t G [r,T], 

where u\ (•) , «2 (•) solutions corresponding to u\ and u 2 , respectively. 

Proof. Arguing as in the proof of Theorem [2] we obtain the inequality 

1 d || \ + ||2 II , n. + ii2 

oil K Ml_M2 ) +a (rti-M 2 ) 



v 



2dt 

< C 3 (t) || (ui - u 2 ) + \\ 2 - I (f 2 (t,v 2 ) - f 2 (t,u 2 ) ,{u\ - u 2 ) + ) dx, 



where v 2 is defined in (1T4T) . 

Using ffT5]) . ^2 < W2) |^2 1 2 < \ui\ 2 + |m 2 | 2 , ([TBI and fl39|) we get 



Thus 

d 
It 



|(wi-W2) + || 2 <2C 3 (t)||(wi-M + || 2 + 2£ / J^(u\-u l 2 ) + dx 



<(2C 3 (t) + l)||(ui-« 2 ) + || +Ke 2 
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for some constant K > 0. By Gronwall's lemma we get 



Mt)-Mt)) + || < 



(u\ - U 2) + || 2 e / r t (2C3( S )+l)d S + Re 2 J' e £(2C 3 (s)+l)ds dr 



<C 2 (r,T)e 5 



Let us consider now the multivalued case. We will obtain first some auxiliary statements. 

We shall define suitable approximations. For any n > 1 we put f^(t, u) = ip n {\u\) f % (t, u)+ 
(1 - ip n (\u\)) g i (t, u), where g i (t, u) = {uf 1 ' u l + f (t, 0, 0), and ip n was defined in ([TBI) . 
Then /„ G C([r, T] x R d ; R d ) and for any A > 0, 



sup sup |/ n (t, m) — f(t,u)\ —7-0, as n — > oo. 

te[T,T] |u|<A 



(40) 



We will check first that f n satisfy conditions fllD-flSD f° r u e where the constants do not 
depend on n. 



Lemma 17 Let f satisfy /or w G R+. For all n > 1 we /iave 

d d d 

53 |£(f,t*)|ifr < ^1(1 + $3 \UT), (fn(t,u),u) >lY,\ U T - D *> 

i=l i=l i=l 

for u G R^ ; where the positive constants D\, D2, 7 do not depend on n. 
if \u\ > n + 1, then for any w G R d we have 

(fnu(t,u)w,w) > 0. 

Moreover, if f,h satisfy (Sty , then f n ,h also satisfies this condition. 

Proof. In view of (El)-© we get 

(fn(t, u),u) = ip n (\u\) (f (t, u),u) + (l- ip n (|u|)) (g (u) , u) 



(41) 



(42) 



u ' 



>i(N) «E ^ - c * + ^ - (N)) E + (* - Vv, (H)) E (*. °> •••> °) 

d d d 

>t(H)«E kr-^+(i - ^ (i«D) 2 E Kr-^i (i - v>„ (i«D) E \f (*> °> •••> 

8=1 

d 



i=l 



i=l 



i=l 
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where a = min{|,a}, for some constant K\ > 0. Also, 

d / d 

E i/to «)i^<ME !/*(*. + E k 



Pi * — , ■ , p i 



pj-1 



i=l \i=l i=l 



< # 3 (c7i(i + e kh + E \ u T + E I/' (*• °> °)K :i 

V i=i t=i t=i 
< k 4 ( E \ U T + 1 

\i=l / 

for some constant K4 > 0. Thus, for D\ = K4, D 2 = C2 + K\C\, a = min{|,a} we have 
( 14T1) for the functions /„. 

Moreover, if \u\ > n + 1, then for any w G M d , 

d 

(fnu(t,u)w,w) = (g u (t,u)w,w) = E (Pi - 1) w ? - °- ( 43 ) 

i=l 

Finally, if ( 157)) is satisfied, then 
W (t, x) - f n (t, u) 

= ^ n (\u\) (h* (t, x) - r (t, u)) + (1 - V n (\u\)) {h* (t, x) - f (t, 0, 0)) > 0, 
for all i, a.e. (t, x) G (r, T) x Q and w such that u l = and > if j 7^ z. ■ 

Let 2 < qi < pi, i = l,...,d. We define also the following approximations f^(t,u) = 
V„(H)/* (t,u) + (l-V n (N))flr*(t,u), where flr*(t,u) = |«f i_2 «H |i*T 1-2 + (t,0,...,0). 
Then (J1D|) holds. We check that /„ satisfy conditions dH)-© f° r u G R+, where the constants 
do not depend on n. 



Lemma 18 Le£ / satisfy /or w G R+. For all n > 1 we have 

d d d 

Ei/to«)i^ <A(i+EKr), (7n(t,«),«) >7Ei^i Pi - D 2, (44) 

8=1 i=l i=l 

/or u G Rj", where the positive constants D\, D2, 7 do not depend on n. 
if \u\ > n + 1, i/ien /or any tu G R d we /tai>e 

fnu{t,u)w,w) > 0. (45) 



Moreover, if f,h satisfy (31), then f n ,h also satisfies this condition. 
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Proof. In view of (141]) we have 

d d 

(f n (t, «),«) = (/ n (t, u), u) + (1 - ^ n (M)) E \u'\ q * > 7 J] \ U T ~ D 2 , 



i=l i=l 



d / d d 

Pi * — ■» . PilH- 1 ) 



E i/to^i 5 * <^i E i/nMi^ + E 

i=l \i=l i=l 



M Pi" 1 



< ^! ^(i + e iw*r) + E i M T i ^ ^ 2 1 + E i M T i ' 

\ i=l i=l J \ i=l ) 

where we have used that pi > qi implies < Finally, ( 14"5)) and condition ( 157)) are 

proved in the same way as in Lemma [T71 ■ 

For every n > 1 consider the sequence (i, w) defined by f^(t, u) = L d p e (s)b n (t, u — s)ds, 

where either b n = f n or b n = f n , defined before. Since any b n are uniformly continuous on 
[t, T] x \—k — 1, k+ 1], for any fc > 1, there exist G (0, 1) such that for all u satisfying 
|w| < k, and for all s for which |« — s\ < we have 

sup \b n (t,u) - b n {t,s)\ < -. 
te{r,T\ fc 

We put /*(*,«) = fn'"(t,u). Then /*(*,-) 6 C 00 ^; for all t G [r,T], fc,n > 1. 
Since for any compact subset Actf and any n we have /* — > 6 n uniformly on [r, T] x A, 
we obtain the existence of a sequence G (0, 1) such that 5 n k — > 0, as k — > oo, and 
(t, -u) — 6^ (t, u)\ < 5 n k, for any i, n and any w satisfying \u\ < n + 2. We define the 
function F% = (F^ 1 , F% d ) given by 

F* i (t,u) = f!i i (t,u)-p8 nk , (46) 

where pGN. 

Lemma 19 Let f satisfy Q)-^ /or w G M+. For all n, k > 1 we have 

d d d 

E i^mi^ < d 3 (i + e m Pi ), (*?(*,«),«) > «/$> 4 i" - ^ (47) 

j=l i=l i=l 

/or m G Ri, where the positive constants D 3 , F> 4 , i/ do not depend neither on n nor k. 
Moreover, if f, h satisfy [Sty, then F*, h also satisfy this condition if \u\ < n + 2. 
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Proof. Since f n satisfy f )4T|) and f n satisfies (jSJ), we have 



d 



i=l \i=l 

d 



< 



^ r i (e ((//^ (s)rfs ) P1_1 Xa (s) \ Unit ' u ~ s) ^ u ds ) + 1 



< R 4 (l + 



for some constant > 0. 
On the other hand, 



ii 

u" 



(F%(t,u),u) = / p €k {s){b n {t,u - s),u- s)ds+ / p 6k {s) (b n (t, u - s), s) ds - p5 nk V] 
^ l/^ s )^it\ ul - S T - AO* - J nd Pe k (s) £ I b\ (t,u- 8 )\£* + R 5 |s'f ) ds 

n „ d „ / d \ 

J ^ d 7=? J ^ d V 2 7=T J 

n „ d n d 

-pSnk ^2 u i > | / p et .( g ) yi k - g T' ^ - #7 -p^yy > ^y 

i=l ^ Rd i=l »=1 »=1 



for some constants u, R% > 0, where in the last inequality we have used that for some D > 0, 
= \ u i - 8 i + gf' < D (\u* - s'^ + js^ 1 ) < D (\u* - sf' + e%) . 



Hence, (@7D holds. 

In view of Lemmas [T71 [T8l the functions b n , h satisfy (1371) . Hence, I/*' (t,u) — b % n (t,u)\ < 
S n k, for any i, n and any u satisfying \u\ < n + 2, implies that 

h* (t, x) - (t, u) = h< it, x) - ft (t, u) + P 6 nk 
>h l {t,x)-b i n {t,u)>0, 

for u such that it, = 0, Uj > 0, j ^ i, and \u\ < n + 2. ■ 
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Define a smooth function (fi n : R + — > [0, 1] satisfying 

( 1, if 0< s<n + 1 + 7, 

0„ (s) = < < n (s) < 1, if n + 1 + 7 < s < n + 2, 
[ 0, if s > n + 2, 

where < 7 < 1 is fixed. Let l*(t, u) be given by 

l*(t, u) = <j> n (\u\) F k n (t, u) + (1 - n (|«|)) b n (t, u) . 

Since for any compact subset AcK d and any n we have /* — > b n uniformly on [r, T] x A 
as k 00 it is clear that 

sup sup \l n (t, u) — b n (t, u)\ — >• 0, as A; —7- 00. 

te[r,T] |m|<A 

Lemma 20 Lei / satisfy /or w G M + . For a// n,k > 1 we have 

d d d 

< + £ 1**1"), K(t, «), u) > A £ |u<|« - D 6 , (48) 

i=l i=l i=l 

for u G R+, where the positive constants D 5 , D 6 , X do not depend neither on n nor k. Also, 

(l k nu {t, u)w, w) > -D?{k, n) \w\ 2 , Vu, w, (49) 

where Dj(k,n) is a non-negative number. 

Moreover, if f,h satisfy (3l\ ), then l k ,h also satisfies this condition. 

Proof. The inequalities given in (j4HJ) are an easy consequence of (|4T|) . (1441) and (1471) . 

On the other hand, if u is such that Ui = 0, Uj > 0, j 7^ i, then in view of Lemmas [TTl 
[TBI and [191 we have 

<j> n (|«|) (tf (t, x) - F^ (t, u)) > 0, (1 - n (|«|)) (tf (t, *) - 6^ (t, u)) > 0, 

as 4> n (\u\) = 0, for \u\ > n + 2. Hence, h l (t,x) — 1% (t,u) > 0, so that condition (1371) holds. 

It is also clear that l k n (t, u) is continuously differentiable with respect to u for any t and 
u. We obtain the existence of D 7 (k, n) such that (jl9~j) holds. Indeed, if |u| < n + 1 + 7, then 
Z*(t,u) = F k (t,u), so that 

|(/L(t, M K W )| = |(Fi(t, M )u;,^)| = \(fH u (t,u)w,w)\ (50) 

/■ 

< \w\ / I Vp e)b (w — s) I |6 n (s, u)\ ds < Rx (k, n) \w\ . 



17 



If \u\ > n + 2, then = b n(t,u), so that by (|42]), (J45]) we have (l* u (t,u)w,w) 

(b nu (t, u)w, w) > 0. Finally, if n + 1 + 7 < \u\ < n + 2, we have 

«K w ) = (N) (^(i, u)w, w) + (l- <\> n (|«|)) (&„„(*, u)w, iw) 



i,j=l 3 



1 



> —R 2 (k, n)\w\ , 

where we have used similar arguments as in (IBTIi) . (H2"j) . (H51) and also that 
9 



(9m 



A(H) 



< i? 3 (n) , w) < i? 4 (n) , F n fci (t, u) < ife (fc, n) 



for any it satisfying n + 1 + 7 < |w| < n + 2, t G [r, T] and any z, j. ■ 

Now we are ready to obtain the weak comparison principle for positive solutions. 

Theorem 21 Let fj,hj satisfy (3o}) -(37). Assume that fj,hj satisfy Uty-^ and (Tl\) for 
u G IR+. We suppose that either f\ or f'2 satisfies $B) for u G and for an arbitrary 
Rq > 0. If < u\ < u\, there exist two solutions Ui,u 2 (of (Tj|) and l[W\). respectively, with 
Ui (r) = u\, u 2 (t) = u^, such that < Ui (t) < u 2 it), for all t G [r, T}. 

Proof. For f~ let us consider the approximations fi n , f 2n defined before with = p\ < pj, 



where the last inequality follows from Lemma [121 Then by Lemmas [TTJ, [18] we have that 
fin, fin satisfy (El)-© for u G with constants not depending on n. Also, (H2l) . ( 1451) hold 
and (J37j) is satisfied in both cases. Moreover, by (JTTT) for u G Mi, we have 



f\ n (t, u) = iJ n (\u\) f[ (t, u) + (1 - iJ n (\u\)) \u 



> v„ (\ u \) r 2 (t, u) + (1 - v„ (M)) (j«T 1 + /2 (*, 0, o)j = r 2n (t, u ) , 

if u G 

As explained before, we can choose a sequence 5 n k G (0, 1) such that 5 n k — > 0, as fc — )■ 



+ kt 



+ • 



+ /i(t,0,...,0) (51) 



00, and 



fin (*> u ) ~ /in (*> M ) < ^nfe, /aJ (*, u) - (*, M ) < ^fc, for any i, n and any u 



satisfying \u\ < n + 2. Further, we consider the functions 



^2n(^«) = /2n(i^)-35n fe . 

By Lemma [19] we know that satisfy (|I])-(jEJ) for u G with constants not depending 
neither on k nor n, and condition (1371) for |u| < n + 2, as well. Moreover, by (1511) we have 

*£(f,u)=/£(f,u)-35„* (52) 
(t,u)-26 nk <fi n (t,«)-2«J nfc 
<f£{t,u)-6 nk = F%(t,u), 
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if u G and \u\ < n + 2. 

Suppose, for example, that f 2 satisfies condition (JHJ) for u G Rl. For any i G [r, T] and 
any m, t> G such that = i>* and u J < f J if j 7^ z, \u\ , \v\ < n, we have 

(«, = An (*, «) - 35nfc > /2n (*, «) ~ ^nk (53) 

= (t, u) - A5 nk > f 2 (t, v) - A5 nk 

= f 2n (t, v) - A5 nk > /« (t, v) - 55 nk = F« (t, v) - 25 nk . 

Hence, ( 13"9~j) is satisfied with R = n and £ = 2<5 nfe . 
Now, we will define the functions 

l k ln (t, u) = n (|«|) F* (t, u) + (1 - n (|«|)) 7 ta (t, it) , 
it) = n (|«|) F 2 k n (t, «) + (1 - n (|«|)) / 2 „ (t, it) . 

By Lemma l20lthese functions satisfy (jlJJ-Q for u G with constants not depending neither 
on k nor n, inequality (Jl9]l . and condition ( 13~T|) . 

In view of (loTj) . (j52|) and n (|it|) = 0, if \u\ > n + 2, we obtain 

n (H) f* (t, «) > n (|«|) f£ (t, u) , (i - n (|«|)) ][ n (t, u) > (i - n (|«|)) A, (t, «) 

and then if u G 

ii»(*.«)>i(*.«)- (54) 

On the other hand, since ^(^ M ) = ^2n (£, w) if |w| < w + 1 + 7, ( 1531) implies 

lL(t,u) >l 2 z n (t,v)-28 nk , (55) 

for any t G [r, T] and any u, v G IR+ such that = t/ and w- 7 < f J if j 7^ z, |tt| , \v\ < n. 
Thus, (1391) is satisfied with R Q = n and e = 25 n fc. 
We consider now the problems 



and 



du 

— -aAu + l* n (t, u) = hj (t, x) , (t, x) G (r, T) x fi, 
u\ X £dn — 0, 



du 

— — a Am + bj n (t, u) = hj (t, x) , (t, x) G (r, T) x O, 



(56) 



— 0, 



(57) 



where 6 ln = f ln and 6 2 n = fin- 

In view of Lemma [201 and [T2| Lemma 5] problem fl56l) has a unique weak solution (•) 
such that u k n (t) > for all t G [r, T]. Let Mj„ (•) be the solutions of ( |56l) corresponding to the 
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initial data u J T , where j = 1,2. Using Lemma [201 it is standard to obtain estimate (fT3l) with 
a constant D not depending neither on n nor k. Then the solutions u^ n (•) of f )56|) satisfy 



it 



?„ (0|| 2 < H\\ 2 + D j (||^(r)f + l)dr = K 2 (||<||, r,T) . 



Thus, since Z£ n satisfy condition fl!| with i? 2 = n 2 > 2 max^ 2 (||^||, r, T) , fT 2 (||w 2 ||, r, T)}, 
by Theorem [161 we know that as u\ < w 2 > we have 



for all t G [r,T], all k and n > (2 max{ir 2 (||^ ||, r, T) , K 2 (\\u% r, T)}) 5 . 



(51 



Arguing as in the proof of Theorem [7] we obtain that for j = 1,2 the sequence vh n 



converges (up to a subsequence) in the sense of (I2"2l - (r2"9~l) to a solution uj n of problem (1BT|) 
with initial data u J T . In particular, as k — > oo we have 



w jn (*) ~> in # for all t G [r, T]. 



(59) 



Fix n > (2max{K 2 (\\u 1 T \\,T,T) } K 2 (\\u 2 T \\ 1 T 1 T)})^, i G and take any t G [r,T\. 

Denote by Q l * k the set 

{xen:v? n {t,x)-v"(t,x)>0}. 



Qit 



By ( )58|) as -> oo we have 

K;(t,x)-<(t,x)) 2 rfx= / ((^(t)-<(t)) + )"rfx<4C 2 (r,T)^,^0. (60) 

Hence, f )59|) implies 

(uin (t, x) — u^n (t, x)) z dx — > I (u\ n (t, x) — u\ n (t, x)Y dx as k — > oc . 



Define the sequence 



A:* 



l ' j I 0,ifxG^V 



Then it is clear that if)\\ L 2^ ~~ >" \\ u in (*) ~~ M 2n (^) 1 1 7- 2(- ^ as A; — )• 00. We note also that 



L 2 {n) 



w n (*) M in (*) — M k (0 weakly in L 2 (Q), as for any £ G L 2 (O) , ([60]) gives 
iu* (t, x) £ (x) 

« (t, x) - < (t, x)) £ (x) dx - [ « (t, x) - < (t, x)) £ (x) cix -> 0. 
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Therefore, w^ 1 (t) — > u\ n (t) — u\ n (t) strongly in L 2 (fl) so that w 1 ^ (t, x) — > u\ n (t, x) — 
u 2n f° r a - a - x e ^- Since w^* (t, x) < 0, for a. a. f2, we obtain 

< ui n (t) < u 2n if) for all t G [r, T\. 

Arguing again as in the proof of Theorem[8]we obtain that the sequences ui n , u 2n converge 
(up to a subsequence) in the sense of (!22l - (!29l) to the solutions Ui,u 2 of problem IQ and 
( ITUl) . respectively. Also, it holds 

< u x it) < u 2 (t) , for all t G [r, T]. 

■ 

As in the previous section we shall generalize this theorem to the case where the constant 
a can be negative. We note that if fj,hj satisfy ( 157)) . then fj (t,v) = e~^fj(t, e^v) + (3v, 
hj = e~P t hj also satisfy (157)) . Arguing as in Theorem [9] we obtain the following. 

Theorem 22 Let fj,hj satisfy ^M)-(3^. Assume that f j7 hj satisfy (EZ|)-(E2P and (122)) for 

u G . Also, suppose that either fx or f'2 satisfies for u G IR+ for an arbitrary Rq > 0. 
If u l < u 1i there exist two solutions ui,U2 (of (G|) and /[TO]) , respectively), with u\ (r) = u\, 
u 2 (r) = u 2 , such that < Ui (t) < u 2 (t), for all t G [r, T]. 



But in this case we can consider another interesting situation, as the values p\ and p 2 are 
not necessarily equal. 

Let /1 satisfy ©-© for mgR{, whereas f 2 satisfy ()50])-(I5D) for !iGl{, with constants 
a\ C{,C 2 , j = 1, 2. Then if a 2 < 0, we make in ([5]) the change of variable v = e _/3i w, where 
P > -a 2 . We obtain problems (152)) -(155)). 

If v (t) is a weak solution of (152)) . then u (t) = e^v (t) is a weak solution of ([9]) (and the 
same is true, of course, for ( 1551) and ( ITUl) ). 

The function /1 (£, t>) = e _/3 */i(t, e^f ) + (3v satisfies (jl])-© for v G with the same p\ 
as /i. Indeed, as e^ 1 ^ 9 ^ 2 ) > 1, we get 

d p\ / d p i d p i \ 

E I /?(*> ^) ^ < ^1 ( E I e N I ^ + E K 1^1 ( 61 ) 

<K 2 (l + J2 



i=l \i=l i=l 

d 

.i\V. 



\v 



1=1 



/ d \ d 

> e" 2/3 ' (a 1 £V' Pi | U T*' +/3|^| 2 > o^^^f -C\. (62) 



8=1 / 8=1 

Then we obtain the following. 
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Theorem 23 Let fj,hj satisfy (f^)-p7j). Let fx satisfy (0)-(G]) for u G R^, whereas f 2 
satisfy [3^) - [3l\l for u G M+, wift constants a^C^C^ j = 1,2. Assume that fj,hj satisfy 
l[Tl\) . Also, suppose that either f\ or f 2 satisfies |2P for u G IR+ for an arbitrary Rq > 0. 
If u t — u 1i there exist two solutions Ui,u 2 (of and [Ttyl . respectively), with u\ (r) = u\, 
u 2 (t) = m^, sitc/j that < Mi (t) < u 2 (t), for all t G [r, T]. 

Proof. We consider problems and (JM]). In view of (lMjl-flgol) and (EU - (JS2D , h (t,v) = 
e- pt fj(t, e^v) + 0v satisfy ®-© for v eR d + with p? = 2, and fx satisfies ®-© for v eR d + 
with the same as fx- Also, defining /ij (t, x) = e'^hjit^x) it is clear that f fTTj) and f[3"T]) 
hold. Finally, if, for example, fx satisfies OH]) for u G for any R > 0, then it is obvious 
that for fx is true as well. 

Hence, by Theorem [2T1 there exist two solutions vi,v 2 (of (132]) and (1331 . respectively), 
with (r) = e'^ul, such that < t>i (t) < t>2 (t), for all t G [r, T]. Thus 

< ui it) = e&vi (t) < e pt v 2 (t) = u 2 (t) , for t G [r, T], 

and ui,u 2 are solutions (of ([9]) and ffTU]) . respectively, such that Uj (r) = u{. ■ 

Remark 24 satisfy ([?]), £/ien i/ie solutions Ux,u 2 given in Theorems\2^ \2^are unique 
for the corresponding initial data. 

Remark 25 All the result proved so far are true if insead of Dirichlet boundary conditions 
we consider Neumann boundary conditions 

^ = in dVt, 
ov 

where v is the unit outward normal. In such a case the space V will be (H 1 (£l)) d . The proofs 
remain the same. 

4 Applications 

We shall apply now the previous results to some model of physical and biological interest. 

4.1 The Lotka-Volterra system 

We also study the Lotka-Volterra system with diffusion 

— = DxAu 1 + u 1 (ax (t) - u 1 - a 12 (t) u 2 - a 13 (t) u 3 ) , 

— = D 2 Au 2 + u 2 (a 2 (t) -u 2 - a 21 (t) u 1 - a 23 (t) u 3 ) , (63) 
du 3 

— = D 3 Au 3 + u 3 (a 3 it) u 3 a 31 (t) u 1 - a 32 (t) u 2 ) , 
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with either Dirichlet or Neumann boundary conditions, where u % = u l (x, t) > and the 
functions (t) , (t) are positive and continuous. Also, Di are positive constants and 
tt C M 3 . The initial data u T belongs to (L 2 (fi)) 3 . 
In this case the functions f\, hi are given by 

I -u l (ai (t) - u 1 - a 12 (t) u 2 - a 13 (t) u 3 ) 
fi (t, u) = -u 2 (a 2 (t) -u 2 - a 21 it) u 1 - a 23 (t) u 3 ) | , h x (t) = 0. 
\-u 3 (a 3 (t) u 3 a 31 (t) u 1 - a 32 (t) u 2 ) 

Uniqueness of the Cauchy problem for this system has been proved only if we consider 
solutions confined in an invariant region (for example, in a parallelepiped T> = {(w 1 , w 2 , -u 3 ) : 
< u l < k 1 } when the parameters do not depend on t) (see [15] and [22]). However, in the 
general case for initial data just in (L 2 (fi)) 3 it is an open problem so far. 

System (IBUj) satisfies conditions dl])-© with p x = p 2 = p 3 = 3 forw G IR 3 [T2J p. 263]. 
Also, it is clear that ([36])-([37D hold. 

We shall compare with the following system 

du ^ 

— = L> 1 Aw 1 + u 1 {a, (t) u}) , 

^ =J D 2 A M 2 + u 2 (a 2 (t)- M 2 ), (64) 

<9li 3 

— =D 3 Au 3 + u 3 (a 3 (t) it 3 ) , 

which is a system of three uncoupled logitic equations. The functions f 2 , h 2 are given by 

f-u 1 (ax (*)-w x )\ 
/ 2 (t, u) = -m 2 (a 2 (t) - u 2 ) , h 2 (*) = 0. 
\-u 3 (a 3 (t)-u 3 )J 

It is easy to see that system ( jMj) satisfies conditions (H])-® with pi = p 2 = p 3 = 3 
fori* G IR 3 . Also, it is clear that (I36 i) -( 13~TI) hold, and that condition §8§ is trivially satisifed. 
Moreover, 

fi(t,u)>fi(jt,u), 

for all t, u G IR+ and i, so that dTTJ holds for u G M 3 . 
Also, we have 

(/2u (*» u ) w, w) > -ai (t) w 2 - a 2 (t) w 2 - a 3 (t) w 3 

+2 (m 1 ^ 2 + m 2 w^ + u 3 wf) > -C \w\ 2 , 

for all uER 3 + ,te [r, T] and w G R 3 , where C > 0. 

Thus, by Theorem (2TJ Lemma [TT] and Remark [25] we obtain the following theorem. 
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Theorem 26 If < u\ < u 2 , there exist a solution u\ of / fff3|) with u\ (r) = u\ such that 
< Ui (t) < u 2 (t), for all t G [r, T] ; where u 2 is the unique solution with u 2 (r) = u 2 of [B~4\ ) 
in the class of solutions satisfying u 2 (t) > for all t. 

This theorem says that there exists at least one solution of the Lotka-Volterra system 
which is dominated by the unique non-negative solution of the uncoupled logistic system 

(El. 

Further, we shall compare with the uncoupled linear system 



dt 
du 2 

du 3 



DiAw 1 + u 1 a 1 (t) , 
D 2 Au 2 + u 2 a 2 (t) , (65) 
D 3 Au 3 + u 3 a 3 (t) . 



Hence, 



-u x ai (t)\ 
f 2 (t, u) = | -u 2 a 2 {t) , h 2 (t) = 0. 
-u 3 a 3 (t)J 

Obviously system (IB3|) satisfies conditions (15Uj) -( 13Tj) forw G M. 3 . Also, it is clear that 
(156"j) -( 1tTT|) and condition ([S]) are trivially satisifed. Moreover, 

fi(t,u)>r 2 (t,u), 

for all t, u G M? + and i, so that (fill holds for u G Also, we have 

{f2u (t, u) w, w) > -a\ (t) w\ - a 2 (t) w\ - a 3 (t) w\ > —C |w| 2 , 

for all ueR 3 ,te [r, T] and w G 1R 3 , where C > 0. 

Thus, by Theorem [23] and Remark [251 we have the following. 

Theorem 27 // < u\ < u 2 ., there exist a solution u\ of I[B1J\) with ui (r) = u\ such that 
< U\ {t) < u 2 {t), for all t G [r, T], where u 2 is the unique solution of (E3]] with u 2 (r) = u 2 . 



Let us consider now the autonomous case, that is, (t) = ai > and r = 0, with Diriclet 
boundary conditions. 

By the changes of variable v 1 (t) = e~ ait u l (t) system (I65p becomes 

% = DiAv\ , = 1,2,3, (66) 
v l \xadn = 0, 
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with initial data Vq = Uq. 

The operator A = -D { A in the space L 2 (Q) with domain D (A) = H 2 (Q) n H] (Q) is 
sectorial |9j. Moreover, since the eigenvalues of A are < -DjAi < Di\ 2 < we have that 
the minimum eigenvalue is strictly greater than 0. Denote by e~ At the analytic semigroup 
generated by the operator A. Then v (t) = e~ At v is the unique solution of ([55]) with 
v (0) = v . 

It is well known [9] that the operator A generate a scale of interpolation spaces X a = 
D (A a ) with the norm \\v\\ a = \\A a v\\ L2 ^, ot > 0, where X a C i^ 2a (fi) with continuous 
embedding. Take < 5 < DiXi. Then by Theorem 1.4.3 in [5] we obtain 

ll^^oL^) < C-r-e"** ||^|| i2(n) for t > 0. 
Since this is true for every 5 < -DjAi we obtain that 

II A a p~ At n i II <n +- a P - D i x it |L,i || 

ll e Atv o\\H 2a (n)— a ^ ° e D * Al< ||' i; o||l2(q) for t > 0. 

We note that H 2a (O) C L°° (O) with continuous embedding if a > |. Since the constants 
C Q are bounded for ck in compact sets, we obtain the existence of C such that 

Ir y o||Loo (n) ^ ^ r 4 e || v o|| L 2(n) 101 r > u > 
where Z) = min{_D!, D 2 , D 3 }. Then the unique solution of (!65|) with it (0) = tt satisfies 

< Ct~'ie {a - DXl)t \\uq\\ fort>0, (67) 



l(L°°(Q)) d 

where a = max{a 1; a 2 , a 3 }. 

Joining ( 167]) and Theorem 1271 we obtain the following result. 

Theorem 28 There exist at least one solution u (t) of the autonomous system KUB\) with 
u (0) = Mo and Dirichlet boundary conditions such that 

< cr'ie {a - DXl)t |K II forte (0,T], 



\\ U \ t )\\(L ao (0)) 3 

where C > and D = min{Di, D 2 , D 3 }, a = max{di, a 2 , a 3 }. 

We shall obtain also a weak maximum principle for the autonomous Lotka-Volterra sys- 
tem with Dirichlet boundary conditions. 

By the maximum principle for the heat equation it is well known (see [I]) that for any 
t > the unique solution of equation ([66]) satisfies 

< v 1 (t, x) < sup Vq for a. a. x G Q. 

n 

Hence, the unique solution of equation ( 165]) satisfies 

< u l (t, x) < e ait sup Vq for e n. (68) 

n 

By ([68]) and Theorem [27] we obtain the following weak maximum principle. 
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Theorem 29 There exist at least one solution u (t) of the autonomous system / Tffgj) with 
u (0) = u and Dirichlet boundary conditions such that for any t e [0,T], 

< u l (t,x) < e ait sup u for a. a. x 6 Q, i — 1,2,3. 
n 

In particular, if Uq G (fi)) 3 , £/ien 

\\ u (*)ll(L°°(fi)) 3 < e a * ||mo|| (loo{q)) 3 /ortG [0,T], 
where a = max{ai, 02, a%}. 

4.2 A model of fractional-order chemical autocatalysis with decay 

Consider the following scalar problem 



r trS +(1 - u) "' 

£(0,t) = £<«,*) = <>, 
aa; ax 

u|t=0 = MO (^) , 



(69) 



where u > 0, N = 1, d = 1, Q = (0, L), and A; > 0, < m, r < 1. The initial data u belongs 
to L 2 (0,a). This equation models an isothermal chemical autocatalysis (see [16]). In [16] 
the authors study the travelling waves of the equation in the case where Q = (0, +00) with 
Neumann boundary conditions at x — 0. The variable u is non-negative, since it represents 
a chemical concentration. 

The funtions /, h are given by / (it) — (u — 1) u m + ku r , h = 0. Clearly, conditions (jl])-© 
hold (with p = m + 2) for it > 0. In this case (I36l) - fl37j) and (jHJ) are trivially satisifed. 

We take /1 = /2 = f, and applying Theorem [21] and Remark [25] obtain the following. 

Theorem 30 If < ul < u 2 ., there exist solutions u±,U2 of $EM) with u\ (0) = u\, m 2 (0) = 
Uq such that < Mi (t) < w 2 (t), for all t G [0, T]. 



4.3 A generalized logistic equation 

Consider the following scalar problem 

' du d 2 u 



_ = _+l 

dt dx 2 

au .„ . au , . 
_(0, i ) = -(a,()=0, 

u|t=o = m (x) , 



(70) 
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where u > 0, N — 1, D, — (0,L), r, q > and r + q > 1. The initial data uo belongs to 
L 2 (0,a). 

This kind of nonlinearities for the logistic equation (instead of the classical (1 — u) u) has 
been considered in [TS1 Chapter 11]. 

The funtions /, h are given by / (u) = {u q — 1) u r , h = 0. Clearly, conditions fll])-© hold 
(with p = r + q + 1) for u > 0. In this case (I3"6]) - (137]) and ([8]) are trivially satisifed. 

We take fi = f2 = f, and applying Theorem [21] and Remark [251 obtain the following. 

Theorem 31 IfO < u\ < u 2 ., there exist two solutionu\,U2 of |70|j withui (0) = Uq, U2 (0) = 
swc/i that < Ui (i) < «2 it), for all t G [0, T]. 
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